Halperin (m, m', n) bilayer quantum Hall states on thin cylinders 
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The Halperin {m,m',n) bilayer quantum Hall states are studied on thin cylinders. In this limit, 
charge density wave patterns emerge that are characteristic of the underlying quantum Hall state. 
The general patterns are worked out from a variant of the plasma analogy. Torus degeneracies are 
recovered, and for some important special cases a connection to well-known spin chain physics is 
made. By including interlayer tunneling, we also work out the critical behavior of a possible phase 
transition between the (331) state and the non-abelian Moore-Read state in the thin cylinder limit. 



A very successful strategy to solve problems in con- 
densed matter physics is to identify the correct simple 
state that a given interacting quantum system can be 
smoothly evolved into. On the other hand, our under- 
standing of fractional quantum Hall states is based on 
radically new techniques Recently, however, a partic- 
ular way has been discussed to adiabatically transform 
fractional quantum Hall liquids into trivial charge density 
wave (CDW) states. In Refs. i, i, i, i it was observed 
that for certain Hamiltonians, a quantum Hall ground 
state does not undergo a phase transition when the two- 
dimensional surface of the system is deformed into a quasi 
one-dimensional (ID) limit, e.g. a thin torus or cylinder. 
In this limit, simple CDW patterns emergeQ. These pat- 
terns are characteristic of the underlying quantum Hall 
state, and can serve as "labels" for the ground states 
and their elementary excitations. Despite their seeming 
simplicity, these labels turn out to be very efficient book- 
keeping tools. They give rise to simple and picturesque 
explanations for some abstract concepts in quantum Hall 
physics, such as the existence of fractional charges and 
topological sectors. Braiding statistics can also be de- 
rived in this language, at least for abelian states ^ll. More- 
over, the efficiency of these ID labels in deriving various 
counting formulas for quasi-hole type states has been sug- 
gested early on by Haldane[8i], who arrived at the same 
effective language by different means d, |^. While such 
formulas can also be obtained by other methods 13, U | , 
Read has recently demonstrated agreement between the 
latter and the ID approachflO]. So far, however, the lan- 
guage of these ID labels has been discussed only for sin- 
gle component quantum Hall states. Here we extend this 
language to the entire class of Halperin (m, to', n) bilayer 
states|12[. To this end, we develop a new technique that 
allows the extraction of the general thin cylinder limits of 
these states directly from the many-body wavefunctions. 
Special attention is paid to the class of (to, to, to) states, 
and to the (331) state. For the latter, implications con- 
cerning the possibility of a continuous phase transition 
into a Moore-Read'13] (Pfafhan) state are discussed. 

Reduction to a ID electrostatic lattice prob- 
lem. In Refs. 0, [^, [^, pseudo-potential Hamiltonians 



have been used to determine the thin torus limits of spe- 
cific quantum Hall states. This method proves to be 
unwieldy for a class as general as the (to,to',7i) states, 
as it would require high orders in perturbation theory 
for large m,m',n. Here we derive a general method to 
extract the limiting CDW patterns directly from many- 
body wavefunctions. Unless otherwise noted, we assume 
TO, to' > n. The {m,m' ,n) states are expected to be 
d — mm'~n? fold degenerate on the torus[14]. Since the 
torus versions of the {m,m',n) wavefunctions arc quite 
complicated, we work on the cylinder, where the same 
CDW patterns must appear between the edges as the 
circumference approaches zero. A straightforward gen- 
eralization of Halperin's {m,m',n) state to cylindrical 
topology leads to the expression 
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vn = (n^l)Q(te};te})xexp(-i 
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where Q(te}; fc}) = U^<Ji^^ ~ iirUi<A^i - 
Cj)™ rij Foi' s = 0, the polynomial part of 

Eq. ([T]) is thus the same that would appear in the usual 
disk topology, with the particle coordinates Za = Xa+iya 
replaced by functions — exp^KZa) periodic in j/q,, 
where k — 27r/Lj, and Ly is the circumference of the 
cylinder. We assume Landau gauge, A = (0,a;), for the 
vector potential, letting the magnetic length equal to one. 
Uppercase (lowercase) indices refer to particles in the up- 
per (lower) layer, whereas Greek indices refer to both lay- 
ers. The integer s in Eq. ([1]) can be interpreted as placing 
s quasiholes in the upper layer at a; — —oo. This does not 
affect the local properties of the incompressible fluid de- 
scribed by Eq. ([1]), which is confined to a finite "ribbon" 
on the infinite cylinder. It is left understood that each 
such ribbon on the cylinder would be infinitely degener- 
ate by translation. We will see, however, that precisely 
d distinct CDW patterns are generated in the interior 
of the ribbon when the thin cylinder limit is taken, and 
when s runs through all possible values and translational 
symmetry is used. 

We will now compute the thin cylinder limit of Eq. ([1]) 
using the following strategy. The lowest Landau level 
(LLL) on the cylinder has a natural basis of ring shaped 
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orbitals V'n = ^"exp(— a;^/2 — K^n^/2) localized at a 
height a; ~ Kn of the cylinder and delocalized around 
the y circumference. We may expand the state Eq. ^ in 
product wavefunctions i/'ni...njv = i^nA^i) ■ ■ ■ • • V'njv(a;w) 
where the N particles occupy definite LLL orbitals. As 
explained in Ref. the coefhcients in this expansion 
are of the form exp(iK^^^n^) • C{ni . . .tin), where 
C(ni . . . n^v) is the coefficient of the monomial ^"^ . . . 
in the polynomial part of Eq. ((T]). From this it is clear 
that those products V'ni...njv i'^ the expansion of Eq. ([T]) 
will dominate as k —>■ oo for which the quantity 

N 
a=l 

is maximal. Since the single particle orbitals ipn are well 
separated along a; as k ^ oo, the dominating products 
'>Pni...nN describe a well defined CDW pattern in the thin 
cylinder limit. To identify this pattern, we observe that 
for any monomial ^"^ . . . that has a non-zero coeffi- 
cient in the polynomial ([liS.!) Q{{£.i}', {£,1}), the expo- 
nents Ua satisfy the following relation: 

na ^^{rriap +Pal3) + s6a.,i (3) 
P 

where niap equals m/2 (m'/2) if both indices are in 
the lower (upper) layer, and n/2 otherwise, Pap is 
an antisymmetric matrix whose values are restricted 
to {— TOq,/3, — mQ,/5 -f 1, . . . , TOQ./3}, and we define Sa./\ = 
<5q,x = J^i^a.i- In Eq. ((S]), each choice for a 
particular value of pap corresponds to the choice of the 
term ^™°'3+p°'3^™q3 Pa/i ^j^g expansion of the factor 

(Ca - ^p?"""^ in Q({6}; {6})- We find that the choice 
of Pais that maximizes Eq. ([2]) is always of the form 

Pa/3 = Sgn(crQ - CT/s) , (4) 

where cr is a permutation of N objects. A proof of this 
statement can be found in Ref. [isj. Figuratively speak- 
ing, Eq. ^ says that each solution of the maximization 
problem Eq. ([2]) corresponds to a "ranking" of the N par- 
ticles by means of the permutation a. If we imagine the 
N particles arranged in a row according to their ranking, 
we may call such an arrangement a state of "the squeezed 
lattice" . The term "squeezed" alludes to the fact that, in 
contrast to the "real space" (orbital) arrangement of the 
particles in the thin cylinder limit that we are seeking, 
there are no empty sites in this squeezed lattice. How- 
ever, we will see below that the particle arrangements on 
the squeezed lattice and in real space are closely related. 
Via Eq. ([4]), the problem is reduced to finding a permu- 
tation (7 that maximizes S. Since particles in the same 
layer are identical, S is invariant under permutations of 
such identical particles on the squeezed lattice. For our 
purposes, the state of the squeezed lattice is thus fully 
described by the data ka = ^s^.t ~ ^&cl^ which equal -1-1 



or —1 depending on whether the particle occupying the 
a-th site, (Tq = (T~^{a), belongs to the upper or lower 
layer. It turns out that in terms of the ka, S can be 
written as 5* = —jn{m + m! — 2n)E, where 

-E = -{P - a)qaqp - 2Qo aqa + const. (5) 

Q</3 a 

and Qa = ka{l + ^qA), A = {m' — m) / {m + m' — 2n) , 
Qo = {2s+N^{ n—m +m—m') I (m-f m'— 2n), 

and we write {N^ ) for the number of particles in the 
lower (upper) layer. Eq. ^ can be interpreted as an 
electrostatic energy assigned to each configuration of the 
squeezed lattice, where a charge q^'^ = ±(1 ± A) is as- 
signed to upper and lower layer particles, respectively. 
The first term in Eq. ([5]) is a linear ID Coulomb inter- 
action between particles. The second term can be inter- 
preted as a linear potential due to an external charge 
distribution, e.g. a charge -|-Qo at the left boundary 
and a charge — Qo at the right boundary of the system. 
Once the charge configuration minimizing E is found, 
the corresponding real space configuration, i.e. the or- 
bital position Ua of the a-th particle, follows from an 
"un-squeezing" rule obtained from Eqs. ([2]), (21) 

^{m'Nl + nNi^s if fc„ - +1 
\mNl^nNl if fc„ = -1, 

where N\ (A^j;) is the number of particles with ka' = +1 
{ka' = —1) to the left of site a on the squeezed lattice, i.e. 
for a' < a. Eq. ([5|) can be viewed as a discrete version 
of the plasma analogy in the thin cylinder limit. 

Solution of the electrostatic problem. We want 
to minimize the energy Eq. ^ for given parameters m, 
to', n, and s. It is useful to write m — n — gr, m' — 
n — gr' , where r and r' are coprime integers. Since 
there are no vacancies on the squeezed lattice, the particle 
configuration and thus the energy E is fully determined 
by specifying the positions of the lower layer particles 
alone. If aj is the position of the j-th lower layer particle 
on the squeezed lattice, we find 

E = q^{q'^ - q'-) ^(a^ - Xj)'^ + const, (7) 
i=i 

r + r' s r 1 

It is now clear how to minimize E and thus maximize S. 
We position the lower layer particles as closely as possi- 
ble to the minima Xj, subject only to the constraint that 
the positions aj are integer, and that 1 < ai < a2 < 
. . . < Qfjvi ^ This determines a pattern of lower and 
upper layer particles in squeezed space. Certain bound- 
ary effects may be present as follows. For example, the 
equilibrium position of the first lower layer particle, xi, 
may be negative. Then for a range of j values, the Xj's 



3 



will be inaccessible by the respective particle coordinates 
Uj . This will lead to a clustering of lower layer particles 
{q^ charges) at the left boundary of the squeezed lattice. 
This in turn can be interpreted as a screening cloud to 
screen the charge +Qo discussed below Eq. ([5]). How- 
ever, for sufficiently large particle number j there will be 
a regime where the particle position aj can be placed 
within a distance 1/2 or less from Xj. This regime of 
the squeezed lattice we will call the "bulk" as opposed 
to the screening cloud on either end of the squeezed lat- 
tice. The appearance of different regimes on the cylin- 
der is expected for the general wavefunction Eq. ([T]) , not 
only in the thin cylinder limit. In particular, one may 
consider Eq. ^ with /N^ much different from the 
correct ratio /v^ of the upper and lower layer filling 
factors in the (m, m', n) phase, where = (rn—n)/d and 
= {m'—n) I d, respectively. In this case Eq. ([T]) is known 
to describe a phase separated state with phase boundaries 
between an (m, m' , n) bilayer phase and a Laughlin-type 
phase in either the upper or lower layer J6J . It is the 
bulk region of the squeezed lattice which corresponds to 
the thin cylinder CDW pattern of the (m, m',n) phase 
when the squeezed lattice is transformed into real space 
via the un-squeezing rules Eq. ([6]). We now derive some 
properties of these patterns. First we note that some of 
the Xj may be half-odd integer. In this case both values 
ckj = Xj ±1/2 lead to a minimum of E. The thin cylinder 
limit of the wavefunction Eq. ([T]) is an equal amplitude 
superposition of the states generated by all these configu- 
rations. In squeezed as well as in real space, the patterns 
corresponding to the two cases aj = Xj H/2 differ by 
an exchange between two adjacent particles in different 
layers. For example, for m = 5, m' — 3 and n — 2 with 
s = 0, we find the following real space thin torus pattern 
from the procedure described above: 

XOXOO T 00 T 00X0X00 t OO t OO . . . (9) 

where characters denote the states of consecutive LLL 
orbitals, meaning unoccupied, and XOX represents a 
triplet t i i t • The unit cell of the pattern Eq. ^ 
has 11 sites. Using translational symmetry, this fully ac- 
counts for the ground state degeneracy of the (532) state 
on a torus. The general situation is slightly more com- 
plicated. From Eq. ([8]), we observe that when j j + r' , 
Xj increases by the integer r -\- r' , and the pattern in 
squeezed space will repeat itself. The squeezed space 
unit cell thus has size r + r' . It consists of r' charges qi 
and r charges q^ and is thus neutral. The size of the unit 
cell in real space can then be obtained from Eq. ([6]) with 
Nl+ r, N^^ r'. It is easy to see that both 
lines in Eq. ([6]) yield Ua ^ na + {mm! — •n?)/g. The real 
space unit cell of the thin cylinder pattern thus has size 
d/g. Note that from this we correctly obtain v"^ = rg/d^ 
i/i- = r'g/d. On a (thin) torus, a given pattern thus only 
accounts for a fraction l/g of the full degeneracy, using 
translational symmetry. However, we see from Eq. ^ 



that varying s will in general yield different real space 
patterns. By "different" , we mean that the patterns are 
not related by translation. The patterns only repeat, up 
to a translation, when s ^ s + g. To see this, we note 
that one can find integers a, b such that ar = r'b -I- 1. If 
we let s ^ s + g, j j + a in Eq. ([5]), we find that Xj is 
shifted by the integer amount a + b. This means that the 
bulk pattern on the squeezed lattice is merely shifted by 
a constant. It is then easy to show from Eq. ((6]) that the 
real space bulk pattern is also just shifted by a constant. 
By varying s, we thus generate g different bulk CDW pat- 
terns in the thin cylinder limit, each having a unit cell of 
size d/g. On the torus these patterns, and those related 
by translation, correspond to g{d/g) = d distinct ground 
states, as expected. We see that ion g > 1, i.e. whenever 
m' — n and m — n are not coprime, not all of the ground 
states on the torus are related by translation. Further- 
more, one may show from Eq. ([8]) that spin fluctuations 
as discussed for Eq. ^ are always present in precisely 
one of the g ground state patterns. In contrast, charge is 
always frozen in the thin cylinder limit, and hence these 
fluctuations are unique to muti-component systems. In 
the following we discuss some cases of special interest. 

(m, m, m) states. For m — m' = n the quantity S, 
Eq. ([2]), does not depend on the permutation cr, Eq. ([4]). 
Every conflguration of upper and lower layer particles 
on the squeezed lattice, for N\ flxed, has the same 
amplitude in the thin cylinder wavefunction. The corre- 
sponding real space configurations obtained from Eq. ([6]) 
live on a diluted lattice where every m-th LLL-orbital 
is occupied. The thin cylinder wavefunction is thus an 
equal amplitude superposition of states on the diluted 
lattice with given N\ NK If the layer index is regarded 
as spin- 1/2 index, these thin cylinder wavefunctions are 
ground states of the ferromagnetic Heisenberg Hamil- 
tonian —J^j{S,nj ■ <S'(m+i)j — 1/4). Using degenerate 
perturbation theory, we have verified that the pseudo- 
potential Hamiltonian whose exact ground state is the 
(to, TO, to) state (see, e.g., RefflTi) does indeed assume 
this form in the thin cylinder limit, to leading order 
in exp(— k^). This shows that the low-energy sector of 
the (to, to, to) pseudo-potential Hamiltonian has a gap- 
less quadratically dispersing branch even in this limit. 
By adding a generic perturbation that breaks the SU (2) 
symmetry down to an easy-plane symmetry, the disper- 
sion of this mode becomes linear, as is well known from 
the theory of spin-1/2 chains. Thus, the low energy 
spectrum we obtain in the thin cylinder limit completely 
agrees with the spectrum expected for (to, to, to) states 
on infinite two-dimensional (2D) surfaces [18j. 

The (331) state. According to the general framework 
established above, the eight (331) ground states on the 
torus come in two classes of four states. These classes are 
distinct by the fact that their respective members evolve 
into one of the following two CDW patterns in the thin 
torus limit, up to translations: 
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iOTOiOTOiOTOiOTOiOTo (lOa) 

XX XX XX XX XX (10b) 

where again XX stands for fj. + |t. We have shown 
that indeed both these states remain zero energy eigen- 
states of the (331) pseudo-potential Hamiltonian (see, 
e.g., Ref. Tt') in the thin torus hmit when leading cor- 
rections in exp(— k^) are considered. As expected (c.f. 
RefsH, S 0, S), it is easy to show that one may form 
domain walls between the two patterns in Eq. (llOp that 
carry the correct fractional charge quantum numbers for 
Laughlin quasi-holes in the (331) state, e.g. -1-1/8 in the 
upper layer and —3/8 in the lower. There is an inti- 
mate connection between the wavefunctions of the (331) 
state and the ~ 1/2 Pfaffian state[l^. This is par- 
ticularly suggestive if we regard the thin torus CDW- 
patterns as the "labels" for the respective ground states, 
as discussed initially. For, if we drop the spin index in 
the (331)-patterns Eq. (fTO|) by simply writing 1 for each 
occupied site, we obtain 10101010 ... for Eq. ()10a|) and 
11001100 . . . for Eq. (fTOb]) . These arejust the labels as- 
sociated with the V = 1/2 PfafRan[5|. Note that the 
latter is only sixfold degenerate on the torus, as the four 
patterns related to Eq. (jlOap by translation collapse into 
mere two patterns when spin indices are dropped. The 
question of the possibility and nature of a continuous 
(phase) transition between the (331) state and the Pfaf- 
fian has generated much interest (20l.l2]|. It is instructive 
to study this phase transition in the thin torus limit. 
To this end, we examine the effect of a small inter layer 
tunneling term Ht added to the pseudo-potential Hamil- 
tonian of the (331) state. We first consider the fate of 
the ground state Eq. (jlOap when Ht is added. The ma- 
trix elements of Ht can be regarded as spin flips. Writ- 
ing H = Hq + Hi + Ht, where Hq contains all terms of 
the pseudo-potential Hamiltonian acting between iden- 
tical or nearest-neighbor LLL-orbitals, Hi all remain- 
ing terms, we can treat Hi + Ht in degenerate pertur- 
bation theory in the thin torus limit. Both states in 
Eq. ([To]) are ground states of Hq, as well as all states 
related to Eq. (llOap by spin flips. In the sector spanned 
by states of the latter kind, we find that to leading or- 
der in exp(— K'^), the effective Hamiltonian is given by 
Heff = Jz E(<^2j<^2j+2 + l)+* Ej <^2j^ ^here the a^f 
are Pauli matrices, and Jz ^ exp(— 2k^). The hopping 
parameter t is assumed to be much less than exp(— k'^/2) 
(the gap of Hq) but may be much larger than J^. H^ff 
is a transverse field Ising model which has a well-known 
phase transition at i = J^. As a result of this transi- 
tion, the two ground states descended from Eq. (|10ap 
and its spin flipped counterpart are replaced by a single 
ground state in the large t phase. The same happens 
in the ground state sector related to Eq. (jlOap by a sin- 
gle site translation. No transition is found in the sector 
related to ground states of the type Eq. (jlObp . It was 
already observed in Ref. [s! that ground states in different 



sectors will in general respond differently to a perturba- 
tion of the pseudo-potential Hamiltonian in the thin torus 
regime. In this regime, the degeneracy of the two types 
of ground states shown in Eq. (|10p is not maintained as 
t increases. Yet if we follow these ground states indi- 
vidually, at large t six of the eight ground states evolve 
exactly into the six thin torus ground states of the 1^—1/2 
Pfaffian, with spins polarized along the x direction. It is 
so far unclear to what extent the behavior described here 
can be extrapolated to the regime of 2D quantum Hall 
phases. However, we argue that any complete theory of a 
phase transition between the (331) state and the Pfaffian 
should reproduce this behavior when the system is de- 
formed into a thin torus. In this sense, an intimate rela- 
tion between this "topological" phase transition and the 
transverse field Ising chain transition may be suspected. 

Conclusion. We have derived a method to obtain the 
limiting charge-density- wave form of Halperin (m, m', n) 
bilayer states on thin cylinders/tori. We found that basic 
properties of these states such as filling factors, torus de- 
generacies and fractional charges can be obtained from a 
one-dimensional variant of the plasma analogy. In some 
important special cases, a connection to the well known 
physics of spin-1/2 chains has been made. This also leads 
to a possible connection between phase transitions in spin 
chains and in topologically ordered states of matter. 
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